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Abstract 



We establish a blowup criterion for the two-dimensional (2D) full compressible 
Navier-Stokes system. The criterion is given in terms of the divergence of the 
velocity field only, and is independent of the temperature. It is almost the same as 
that of barotropic case. 
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1 Introduction 

The motion of compressible viscous, heat-conductive, ideal polytropic fluid is gov- 
erned by the full compressible Navier-Stokes system. Suppose that the domain occupied 
by the fluid is f2. The whole system on (0, T) x Q, consists of the following equations 



In this paper, we consider the two-dimensional case, i.e., f2 is a bounded smooth domain 
in IR 2 . p,u = (u 1 , u 2 ) , 9 and P = Rp9 (R > 0) represent respectively the fluid density, 
velocity, absolute temperature and pressure. In addition, D(u) is the deformation 
tensor: 



p t + div(pu) = 0, 

(pu) t + div(pu ®u) — pAu — (p + A)Vdivu + VP = 0, 
c v [{p9) t + dw(pu9)] - kA9 + Pdivu = 2p\®(u)\ 2 + A(divu) 2 . 




together with the initial-boundary conditions 




(1.2) 
(1.3) 




The constant viscosity coefficients p and A satisfy the physical restrictions: 



p > 0, p + A > 0. 



(1.4) 
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Positive constants c v and n are respectively the heat capacity and the ratio of the heat 
conductivity coefficient over the heat capacity. 

There is a huge amout of literature on the existence and large time behavior of 
solutions to (jl.ip . For the case that the initial density is far away from vacuum, local 
existence and uniqueness of classical solutions were proved in [2ll[26][28] • Matsumura- 
Nishida [25J first obtained the global classical solution when the initial data is close 
to a non- vacuum equilibrium in some Sobolev space H s . Later, Hoff [8] constructed a 
global weak solution for discontinuous initial data, with the assumption that the initial 
density is strictly positive. 

Normally the presence of vacuum state makes the problem more complicated. We 
recall some results about the weak solution for this case first. For the global weak 
solution to the barotropic case, the major breakthrough is due to Lions [24] and subse- 
quently improved by Jiang-Zhang [20] and Feireisl [7]. They succeeded in constructing 
a global weak solutionwith finite energy when the pressure P = p 7 , 7 > y, where N 
is the dimension. Recently, Huang-Li [12] obtained the global weak solution to the full 
Navier-Stokes system (jl.ip provided the initial energy is suitably small. Strong solu- 
tions have also been under investigation. The first local existence result was derived 
by Cho-Kim [4J. For a global classical solution with small energy, refer to |12jll6j. On 
the other hand, Xin [33] first contributed a remarkable blow-up result. He proved that 
if the initial density has compact support, any smooth solution to the Cauchy problem 
of the full Navier-Stokes equations without heat conduction blows up in finite time. In 
this direction, for more recent progress, sec [3,27J and references therein. 

Taking into consideration both the local existence results and the blow-up results, 
then it is important to study the mechanism of possible blowup and structure fo possible 
singularity. For the blowup criterion for the compressible flow, there have been many 
results, |5|51BlinfnfI5HI5IIgf^ It should been mentiioned here that Huang-Li- 

Xin [15] first established a Beale-Kato-Majda type blowup criterion for the baratropic 
case. In fact they proved that if T* is the maximal time of existence for local strong 
solution, then 



lim 



I \\Vu\\L°°dt = 00, (1-5) 
* Jo 



under the assumption 7/j, > A when Q is a three-dimensional domain. Jiang-Ou [19] 
extended this criterion to the full Navier-Stokes system (jl.ip over a periodic domain or 
unit square domain of M 2 and proved that 

f T 

lim / \\L°°dt = 00. (1.6) 

t^t*J 

Recently, Huang-Li- Wang [13] obtained a Serrin type blow up criterion for (jl.ip in R^. 
Here is the criterion, 

fT 



lim 



r 2 n 

/ ( \\diru koo + \\u \\ s rrdt) = 00, — I = 1, iV < r < 00. (1.7) 

Jo s r 

which is analogue to the Serrin criterion for the 3D incompressible Navier-Stokes equa- 
tions. In particular, for N = 2, if one can bound a priorily ||k|| j / 2(o i T;L°°)" norm or 
IMlL 4 (o,T;L 4 )" norm > then (jl.7p can be replaced by 

f T 

lim / I diva I roodt = 00. (1.8) 
t^t*J 
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If (]1.8p is proved, it is an improvement of the work by Jiang-Ou [19] and it reveals some 
connection between the compressible and incompressible Navier-Stokes equations, since 
global strong solution has been proved for 2D incompressible case. The question is we 
can not get the uniform bound of |MIl 2 (o,T;L°°) or ll n lli 4 (o,T;L 4 ) from the a priori energy 
estimate. The aim of our paper is to verify (|1.8j) and the key trick is the use of Lemma 
12.21 below, one critical Sobolev embedding inequality. 

The results such as (jl.5p or (jl.6p or (II. 7p , notice that they are all in terms of velocity 
field only. There is another big class of results which are in terms of density p and 
temprature 9. For example, Sun-Wang-Zhang [30J obtained the following criterion in 
3D, 

lim sup {||p||l°° + + ||0||l°°} = oo. (1.9) 

T->T* \o<t<T J 

Fang-Zi-Zhang [6j extended the result to the 2D problem with a refiner form, 

lim sup {||p||l°° + = co- (1-10) 

T^T* <i<T 

Before stating our main result, we first explain the notations and conventions used 
throughout this paper. We denote 

/ fdx = / fdx. 
J Jn 

For 1 < p < oo and integer k > 0, the standard homogeneous and inhomogeneous 
Sobolev spaces are denoted by: 

LP = LP(Q), W k 'P = W k 'P(n), H k = W k > 2 , 
Wq' p = ue W x ^\ « = 0on dQ, Hi = Wq' 2 . 

Let 

f:=f t + u-Vf 

denote the material derivative of /. 

Since we are going to work with the blowup criterion of the strong solutions, we'd 
like to recall the result for the existence of the local strong solution. The solution to the 
3D full Navier-Stokes system with vacuum was obtained by Cho-Kim [3]. The method 
there can be applied to the case in this paper, i.e. the case that Q is a bounded domain 
in R 2 . And the corresponding result can be stated as follows(or refer to [6]): 

Proposition 1.1 Let q 6 (2,oo)&e a fixed constant. Assume that the initial data satisfy 

po>o, poeiy 1 ' 9 , uoe^ntf 2 , 0q s h 2 , 

with the compatibility conditions 

pAu + (A + ^)Vdivn - RV{p 9 ) = pf Si, (l.H) 
kA8 + ^\Vu + (Vn )' r | 2 + A(divno) 2 = pf 5 2, (1-12) 

for some g\,g2 £ L 2 . Then there exist a positive constant Tq and a unique strong 
solution (p,u,8) to the system &l.l\) - (T73\) such that 

p>0, pGC([0,T ];W^), 6 eC([0,T ];H 2 ), (1.13) 
ueC([0,T ];HlnH 2 ), (u, 6) G L 2 ([0, T); W 2 ' q ), (1.14) 
(ut,0 t ) e L 2 {[0,To];H x ), {y/m, y/pOt) 6 L°°([0, T ); L 2 ). (1.15) 
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For the blowup criterion for the local strong solution, here is our main theorem. 

Theorem 1.2 Suppose the assumptions in Proposition \1.1\ are satisfied and (p,u,9) is 
the strong solution. Let T* be the maximal time of existence for that strong solution. 
If T* < co, then 

\\m f ||divu|| L i (0iT;LO o) = oo. (1.16) 

A few remarks are in order, 

Remark 1.1 It is worth noting that the conclusion in Theorem \1.2\ is somewhat sur- 
prising since the criterion (j!.16j) is independent of the temperature and is the same 
as that of barotropic case( [14tf ). In fact, it seems that the nonlinearity of the highly 
nonlinear terms \D(u)\ 2 and (divu) 2 in the temperature equation is stronger than that 
of div(pu (g) u) in the momentum equations ( f31f). however, (|1.16p shows that the 
nonlinear term \Vu\ 2 can be controlled provided one can control div(pu ® u). 

Remark 1.2 It is well known that the 2D incompressible homogenenous Navier-Stokes 
system has a unique global strong solution if the initial velocity belongs to L 2 or some 
more regular space, and recently it is proved in [11] that the 2D incompressible non- 
homogenous Navier-Stokes system also has a unique global strong solution under some 
compatibility conditions, so the result in our paper is reasonable from this point. The 
blowup criterion here shows that dfvu plays an important role in the fluid dynamics. 

The rest of the paper is organized as follows: In Section 2, we collect some elementary 
facts and inequalities. The main result, Theorem II. 2 \ will be proved in Section 3. 



2 Preliminaries 

In this section, we recall some known facts and elementary inequalities that will be 
used later. 

The first proposition is for the Lame system, which comes from the momentum 
equation (ILip ?. Assume that f! C I 2 is a bounded smooth domain. Suppose U S Hq 
is a weak solution to the Lame system, 

liAU + (p + A)Vdiv*7 = F, in ft, 
U(x) = 0, on dn. 

The system is an elliptic system under the assumption (II. 4D . hence some regularity 
estimates can be derived. For a proof, refer to [30] . 

Proposition 2.1 Let q G (l,oo). Then there exists some constant C depending only 
on \,fi,p and SI such that 

• ifFe L p , then 

\\U\\ W 2, P < C\\F\\ LP ; (2.2) 
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ifF G W~ x * (i.e., F = div/ with f = (^)ax2, /y G tten 

ll^llwna. < CH/Hlp. (2.3) 

Moreover, for the endpoint case, if G L°°nL 2 , then VJ7 G BMO(Q) and there 
exists some constant C depending only on A, fi, suc/i t/iai 

\\VU\\ BMO <n)<C{\\f\\L~ + \\f\\i?). (2-4) 
#ere ||g|| B MO(n) : = llfflU 2 + [s]bmo(Q)i 

^Wofn) : = sup 1 / -ffn r (a;)|dy, 

xGH, re(0,d) l s< rWI J!! r (i) 



f2 r (:r) 



where Q, T {x) = B r (x) n f2 and |O r (x)| denotes the Lebesgue measure of £l r (x). 

Two logarithmic Sobolev inequalities will be presented, which originated from the 
work due to Brezis-Gallouet pQ and Brezis-Wainger [2]. The first one, together with 
Proposition 12.11 will give the estimate of Vp. For its proof, see also [30J. 

Lemma 2.2 Let Q be a bounded Lipschitz domain inW? and f G W 1 ' p withp G (2,oo), 
there exists a constant C depending only on p such that 



< 



C(l + ll/l| B MO(n)ln(e+||/|| wl , P )). (2.5) 



The second inequality is in terms of both space integral and time integral. The proof 
can be found in [17J or refer to [23] for a similar proof. It plays an important role in 
the proof of Lemma 13.31 



Lemma 2.3 Let $7 be a smooth domain in R 2 . and f G L 2 (0,T;Hq n W 1,p ), with 
p > 2 . Then there exists a constant C depending only on p such that 

L2(0,T;L~) < C (l + II/IIl2(0,T;/P) ln ( e + II / II L 2 (0,T;W^(Q) )) • (2-6) 



3 Proof of Theorem 11.2 



Let (p, u, 9) be a strong solution described in Theorem 11.21 Suppose that (|1.16p is 
false, i.e., 

lim [ | divit ] | £,i (o,T; ) ^ Mo < +oo, (3.1) 

which together with (11. ll^ yields immediately the following upper bound of the den- 
sity(see |14t Lemma 3.4]). 

Lemma 3.1 Assume that (|3.ip holds. Then it is true that for < T < T*, 

sup HpIIl 00 < C, (3-2) 

0<t<T 

where and in what follows, C, C\, C2 and C3 denote generic constants depending only 
on Mq, ji, A, R, k, c v , T* , and the initial data. 



The next estimate is similar to an energy estimate. 
Lemma 3.2 Under the assumption f)3. 1 1) . it holds that for < T < T* , 



sup [ (p9 + p\u\ 2 ) dx + [ \\Vu\\ 2 L2 dt < C. (3.3) 

0<t<T J Jo 

Proof. Applying standard maximum principle to (jl.lP o together with 9q > (c.f. 
[?,?]) shows that 

inf 9(x,t)>0. (3.4) 

R3 X [0,T] 

It follows from (jl.ip that the specific energy E = c v 8 + ^\u\ 2 satisfies 

(pE) t + d\v{pEu + Pu) = A(k6 + ^p\u\ 2 ) + ^div(u • Vu) + Adiv(udivu). (3.5) 
Integrating (|3.5p over 0, x [0, T] directly yields 

sup (pO + p\u\ 2 ) dx < C. (3.6) 

0<t<T J 

Multiplying (jl.ip o by u and integrating the resulting equation over we obtain after 
using ([33D and (jBU]) that 

1 d f 

-— / p\u\ 2 dx + //||Vu||^ 2 + (/i + A)||divu[|^ 2 
<C||divu||ioo / p^dx 



<C||divu| 



which together with (13. ip and (I3.6P gives (I3.3P . It completes the proof of Lemma 13.21 

For a slightly higher order estimate, we derive the bound for the L°°(0, T; L 2 )-norm 
of Vu, which will play a key role in obtaining higher order estimates. 



(3.7) 



Lemma 3.3 Under the condition (|3.ip . it holds that for < T < T* , 

sup [ (p6 2 + \Vu\ 2 )dx+ [ [ (\V9\ 2 + e\Vu\ 2 + p\u\ 2 )dxdt<C. 

0<t<T J Jo J 



Before the proof of Lemma 13.31 we present an auxiliary lemma, which controls LP- 
norm of 9 by ||V0||^2. And it will help the proof of Lemma 13.31 

Lemma 3.4 Under the condition (|3.ip . it holds on [0, T*) that for every p € [l,oo), 

\\e\\ LP <c + c( P )\\ve\\ L 2. (3.8) 

Proof Denote by 9 = pr f 9dx, the average of 9, 



\9\ J pdx < 



p9dx 



+ 



p(9-9) dx 



(3.9) 



<C + C||V0|| L2 , 
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which together with Poincare's inequality implies 

||0|| L2 <C + C||V0|| i2 . 

And consequently, (|3.8j) holds. 



(3.10) 



Proof of Lemma 13.31 First, multiplying (jl.ljU by 9 and integrating the resulting 
equation over Q lead to 



d_ 

dt 



J P 9 2 dx + 2k\\V9\\ 2 l2 < C\\divu\\ La ° J P 9 2 dx + C j 9\Vu\ 2 dx. (3.11) 



To make the estimate close, one needs to bound the term J 9\S7u\ 2 dx in (|3.1ip . To 
achieve that, we borrowed the idea from [13], by multiplying (11 . 1 by u9 and integrating 
the resulting equation over SI is done. Then 

fj, J \Vu\ 2 9dx + (fj, + A) J \divu\ 2 0dx 

= — J pii- u9dx — p, j f u ■ Vn • V9dx — (p + A) J divu • (u ■ V9)dx — J VP • u9dx. 

(3.12) 

We estimate the terms on the righthand of (I3.12p , By the Young's inequality, 



pu ■ u9dx 



< r] J p\u\ 2 dx + Cr) J , 



q2i ,|2 



u\ dx, 



and 



p j u-Vu- V9dx + (p + A) J divu ■ (u ■ V9)dx 



<7l|V0||i2 +C t j \u\ 2 \Vu\ 2 dx, 



(3.13) 



(3.14) 



where r], e are small positive constants to be determined later. Using integration by 
parts, 



\7Pu9dx 



J P9divudx + j Pu ■ VVd.r 



e 

4' 



< tI|V6I||| 2 + C||divn|| L o 



(3.15) 



dx + C € I p 2 9 2 \u\ 2 dx. 



Combining the estimates (|3.1ip - (|3.15p . we obtain after choosing e suitably small that 
> 2 dx+ I (9\Vu\ 2 + k\V9\ 2 ) dx 



d 
dt 



< Cr] J p\u\ 2 dx + C\\divu\\ L °° J p9 2 dx + C v J (p9 2 \u\ 2 + \u\ 2 \Vu\ 2 ) dx. (3.16) 

Note that there are some terms such as J p\u\ 2 dx whose estimates are not clear. These 
terms look less frightening than f 9\X7u\ 2 dx, if one is familiar with the techniques for 
regularity estimates of compressible Navier-Stokes equation. One standard way is to 
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multiply (|l.ip 9 by u t and integrate the resulting equation over $7. Then by the Young's 
inequality, we obtain that 

\~dt\ ^' Vu ' 2 + ^ + A )( divn ) 2 ] dx + / Pl^dx 
= j P u. { u.V)udx + jp^u t dx 

- \ J P\u\ 2 dx + C J \u\ 2 \Vu\ 2 dx + ^ J Pdivudx- J P t divudx (3.17) 

To deal with the term j Ptdivudx, we employ some technique which is a combination 
of those from [30] and [32]. First, we split u into two parts, v and w. Let v = C~ l VP 
and w = u — v. As noted in |30| . divu; acts as the effective viscous flux for the bounded 
domain case. Now 

J Ptdivudx = J Ptdivvdx + J P t divwdx. 

Herein 



Ptdivwdx = — I VPtvdx 

(Cv)tvdx (3.18) 



and according to (|3.5 
Ptdivwdx 



R 

(. i j 

R 

c 



(pE)tdivwdx — j -(p\u\ 2 )tdivwdx 



— If (pEu + Pu — kV9 — pVu • u — pu ■ Vu — Audivu) • Vdivwdx 

V [J 

+-div(pu)|u| 2 divu)(ix — J pu t udivwdx^ 
=— | J \{c v + R)p6u + -p\u\ 2 u — kV6 — pVu ■ u — pu ■ Vu — Audivu^ • Vdivwdx 

— J -pu\u\ 2 V div wdx — J pu ■ udivwdx^ 

= — { f \{c v + R)pOu — kV0 — uVu ■ u — pu ■ Vu — Audivu] • Vdivwdx 
c v [J 

— J pu ■ udivwdx^ . 

(3.19) 

By virtue of Proposition 2.1, we have 

HVwIIra < C\\ P e\\ L 2, (3.20) 

and 

\\D 2 w\\ L 2 < C\\pu\\ L 2. (3.21) 



Making use of these two inequalities (|3.20p and (I3.2ip . we obtain that 

- J Ptdivwdx < C {\\y/p0\\ L 2\\u\\ LO o + \\V8\\ L2 + \\Vu\\ L 4u\\ L ™) ||Vdivu;|| L2 

+ C||- v /p'u|| i 2||'u||i,oo||divi(;|| L 2 

< CiUv^ll^llullioc +Cs\\V9\\ 2 L2 + C 5 \\Vu\\ 2 L2 \\u\\ 2 Laa + 6\\^pu\\ 2 L2 . 

(3.22) 

Substituting (|3.18p and (|3.22p into (|3,17p . we obtain after choosing 6 suitably small 
that 



d_ 
dt 



j (^i\Vu\ 2 + (p + A)(div-u) 2 - 2Pdi™ + ^ + - + 1 ^ ((-L) 1 / 2 ^ 2 ^ (fx + ^ p|u| 2 dx 
< Ci«||V0|||a + Ci J 9\Vu\ 2 dx + C J p 2 \u\ 2 6 2 dx + C J \u\ 2 \Vu\ 2 dx 

+C (||divn|| L oo + 1) ^|| Vn|| 2 2 + J p6 2 dx^J + C. (3.23) 

Choose some constant C2 > C\ + 1 suitably large such that 



/i|Vu| 2 - 2Pdivn + C 2 p9 2 > ^|Vu| 2 + P 9 2 , (3.24) 
adding (|3. 16[) multiplied by C2 to (|3.23p . we have after choosing rj suitably small that 

^ f ( /i|Vn| 2 + ( M + A)(divu) 2 - 2Pdivn + — ^— + C 2 ^ 2 + h(-£.)V* v A dx 
at J \ 2/z + A 2 / 

+^ / o|u| 2 (fa; + y #|Vu| 2 c&; + k J \V6\ 2 dx 

< C(||divu|U« + 1 ) fl|Vu||£a + Jp6 2 dx^j+C J ( P 6 2 \u\ 2 + \u\ 2 \Vu\ 2 ) dx + C 

< C (||divn|| L oo + 1) • (j|Vu|| 2 2 + j P 6 2 dx^j + C 

+C7||u||ioo (J p0 2 dx + J \Vu\ 2 dx^j . (3.25) 

Let 

¥(t) = e+ sup f ||Vn(r)|| 2 2 + I p9 2 { T )dx\ + f f (p\ii\ 2 + 6\Vu\ 2 + \V6\ 2 ) dxdr. 
re[0,i] V J J JO J 

(3.26) 

By virtue of Gownwall's inequality, for every < s < T < T* , 

9(T) < Ctf(a) explcj^ || u(r) || drj . (3.27) 

Now it is time to get a good control of ||i*||z, 2 (s,T , ;i, o °) - Making use of Lemma 2.3, we 
can get that 

Wh(s,T;L°°) < C (l + \\ u \\h(s,T;Hl) ln ( e + II n lli 2 (s,T;W^)) ) • (3-28) 



(I 
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By Proposition 2.1 and Lemma [3 



\\u 



(3.30) 



1^1,3 <C||p«|| 6 + C\\P\\ L 3+C\\u\\ L 2 

(3 29) 

< C\\pu\\ L 2 + C\\V9\\ L 2 + C\\Vu\\ L 2 + c 

which implies that 

IMIl^TW 1 - 3 ) < C||/5^IIl 2 ( S ,T;L 2 ) + C\\V9\\ L 2( S>T . L 2} + C\\Vu\\ L l( StT . >L 2) + C 

< C^iT) 
Substituting (f3T30l) to (ET2B1) . 

IMIW;^) < C (l + h||i 2{S)T;Hl) In (C*(T))) . (3.31) 

Taking this inequality (|3.31|) back to (|3.2T[) . then we get 

\P(T) < C^{s){C^{T)f iM ^ 2 ^H^ (3.32) 

Recall the energy like estimate (|3.3p . and choose some s which is close enough to T* 
such that 

„ „2 1 

T ig?._ C 3||«|| La(8iT . ff i < 

then 

(T) < C*(s) 2 < oo. (3.33) 

The following Lemmas l3.5M3.10l will deal with the higher order estimates of the so- 
lutions which are needed to guarantee the extension of local strong solution to be a 
global one under the conditions (|3.ip and (|3.7p . 



Lemma 3.5 Under the condition (J3JJ), it holds that for <T < T* , 



o<t<T . 



CI 


(p9 2 + \Vu\ 2 ) 







Proof. First, applying u J [<9 t + div(u-)] to (jl.lP o and integrating the resulting equation 
over Q, we obtain after integration by parts that 

\jt I p ^ dx = ~ f iij [d j Pt + dw(ud j P)]dx + fi J u j [Au{ + div(uAu j )]dx 

3 
i=l 

We get after integration by parts and using the equation (jl.ip that 
N\ = — J u j [djP t + div(djPu)}dx 

= «/ - . w - v, - e^u) * + / 

= R J djV? (^p6 — Opdivuj dx + J Pdivudivudx — J Pd^v? dju k dx (3.36) 

< |||v«||^, + c\\ P e\\ 2 L2 + c J P 2 e 2 \vu\ 2 dx 

< ^l|Vn|| 2 2 + C\\p9\\ 2 L2 + C\\6\\it + C\\Vu\\i 4 . 
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Integration by parts leads to 

N 2 = n J u j [Au{ +div(uAu j )]dx 



p / ( diii 3 diu\ + Av?u ■ Vn J ) dx 



= -A* / (l v ^! 2 - diU j u k d k d iU j - diu j diu k d k u j + Au j u ■ Vn J ) dx (3.37) 

<-y y |Vu| 2 <ix + C /" |Vu| 4 dx. 
Similarly, we have 

iV 3 < -g(A* + A)||divu|| 2 2 + |Vu| 4 cte. (3.38) 

Substituting (f3736]) - (f3738] ) into (13351) implies 

j t J p\ii\ 2 dx + fi\\Vu\\ 2 L2 < C J P 6 2 dx + C\\e\\i 4 + C j \Vu\ A dx + C 

<C J P 9 2 dx + C\\V0\\l 2 + C\\^pu\\\ 2 + C, 

where we have used Proposition 2.1 and Lemma l3.4[ 

Next, multiplying (|l.l|h by 9 and integrating the resulting equation over yield 
that 

c v J p\0\ 2 dx + | (|| V0||| 2 ) t = -k Jv9- V(u • V0)dx + A J (divu) 2 9dx 

+ 2p J \Q{u)\ 2 6dx - R j p9divu6dx ^ ^ 



(3.39) 



4 

A 



i=l 



We estimate each = 1, • • • ,4) as follows: 

First, it follows from Sobolev embedding theory that for any e G (0, 1], 



2 |Vu| 2 dx < C||0||£,»||Vu| 



L 2 



<e||V 2 0|| 2 :2 + C £ ||V0||i 2 +C7 5 , 



(3.41) 



which together with the standard W /2 ' 2 -estimate of gives 

||^|||-2 <C j P 9 2 dx + C J p 2 6 2 \Vu\ 2 dx + C j |Vu| 4 cfx + C||0|| 2 2 

<c J pe 2 dx + Ce\\v 2 e\\ 2 L2 + c e \\ve\\ 2 L2 + c J \v u \ 4 dx + c e . 

Hence, choosing some e small enough, we have 

\\0f H 2 <C J p9 2 dx + C|| W|| 2 2 + C j \Vu\ A dx + C. (3.42) 
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(3.43) 



Consequently, 

|iil<c J \vu\\ve\ 2 dx + c j \ve\\v 2 e\\u\dx 

^cuvuii^iivflii^ + ciiv^ii^iiveii^iiuiu* 
<5\\v 2 e\\ 2 L2 + c 5 \\ve\\ 2 L2 

<Cdj p6 2 dx + C 5 \\V8\\ 2 L2 +Cs j \Vu\ 4 dx + C 5 . 
Next, integration by parts yields that, for any 77, S € (0, 1], 
J 2 =A /" (divu) 2 ^dx + A /" (divu) 2 u • V0dx 

=A ^ (divu) 2 6><ix^ - 2A J 6>divudiv(ti - u ■ Vu)dx + A J (dwu) 2 u ■ V9dx 
=A ^y (divu) 2 6dx^j — 2A J Odivudivudx 

+ 2A y 0divudiv(ii • Vu)dx + A j (divufu ■ VOdx 
=A ^y (divu) 2 6dx^j — 2X J Odivudivudx 

+ 2A y BdivudiVPdju'dx + X j u-V (fl(divu) 2 ) dx 
<A (^J {divufedx^j + C||0|Vu| || L 2 (II Vu\\ L z + I) Vu|| 2 4 ) 

<A ^ (divu) 2 0cfa^ + 77||V«||| 2 + C5 j P 6 2 dx + CfcJ W|| 2 2 + Q||Vu||^ + C s , v , 

(3.44) 

where in the last inequality, we have used ()3.41 j) and (|3.42|) . 
Then, similar to (|3.44p . we have that, for any 77, «5 G (0, 1], 

h<2fi(J \®{u)\ 2 0dx\ + rj\\Viif L 2 + CS J P 9 2 dx 

+ c s jvef L2 + c s \\v u \\% + c s . 

Finally, it follows from ([HIT]) and (|3~4^) that 

\h\ <5 j p9 2 dx + C s y 9 2 \Vu\ 2 dx 

<C5 j p9 2 dx + C s \\Ve\\ 2 L2 +C s y \Vu\ A dx + C s . 



(3.45) 



(3.46) 



Substituting (|3.43p - (|3.46p into ()3.40p . we obtain after choosing 5 suitably small that, 
for any 77 G (0, 1], 



I y (l\ve\ 2 -e[\(divu) 2 + 2pMu)\ 2 ])dx + ^ J P e 2 dx 

< Crj\\Vu\\ 2 L2 + C„||Vu||i 4 + c v \\ve\\ 2 L2 + c v 

< C V \\Vu\\ 2 L2 + Cjp^uWt, + CjV0||| 2 + C v . (3.47) 
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Noticing that 

' [X(divu) 2 + 2fi\D(u)\ 2 ] dx 



<C\\9\\ L 4Vu tt2 



L 12/5 

4/3,1-, , 1.2/3 



<C (\\Ve\\ L 2 + 1) • \\Vu\\% ||Vn||2 4 (3-48) 

<c (\\ve\\ L , + 1) (\\Vpu\\% 3 + ||p0[| J? + 1' 



so adding (|3.39p multiplied by 2r/ 1//2 to (|3.47p , we obtain (|3.34p after choosing r\ suitably 
small and using Gronwall's inequality. Thus we complete the proof of Lemma 13.51 

As a corollary, we can bound ||0||x,4 and ||Vti||^4. 
Corollary 3.6 Under the condition (|3.1|) . it holds that for < T < T* , 

sup (||0|| i4 + ||Vu|U*) < C. (3.49) 

0<t<T 

Proof. By virtue of Lemma 13.41 and Lemma 13.51 

\\9\\ L i <C\\V9\\ L 2 + C <C. (3.50) 

Consequently, according to Proposition 2.1, Lemma 13.41 and Lemma 13.51 

\\Vu\\ L 4 < C\\pu\\ L 2 + C\\ P 9\\ L 4 + C\\Vu\\ L 2 +C <C. (3.51) 

Next, we will derive the desired estimates for 6. In fact, we have 



Lemma 3.7 Under the condition (|3.ip . it holds that for < T < T* , 



sup I p6 2 dx+ [ \\V9\\ 2 L2 dt < C. (3.52) 

0<t<T J Jo 

Proof. Applying the operator dt + div(u-) to (jl.ip a. we use (Jl.ip i to get 

CvP (d t e + u ■ w) 

= kA0 + k (divuA6 - d t {dm ■ V8) - d. L u ■ Vdtf) 
+ (A(divn) 2 + 2^|D(u)| 2 ) divu + Rp6d k u l diu k (3.53) 
— RpOdivu — RpBdxvu + 2A ( divu — d k u l diu k ) divu 



+ p(diu 3 + dju 1 ) [diu 3 + dju 1 - diU k d k u 3 - djU k d k i 
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Multiplying (I3.53D by 9, we obtain after integration by parts and Corollary 13.61 that 
' ( f p\9\ 2 dx) +k\\V9\\ 2 l2 



2 



<C J \Vu\ (\V 2 9\\9\ + |V0||V0|) dx + C j \Vu\ 2 \9\ (\Vu\ +6)dx 
+ C [ p\9\ 2 \Vu\dx + C [ p9\Vu\\9\dx + C [ \Vu\\Vu\\9\dx 



<C||V U || L 4||V 2 0|| i2 ||^|| L 4 + C||Vn||| 4 ||0|| L 4 

+ C||V«||i4||6»||x4||fl|U*+C'||V«|U4|| V ^|| L 2||e|U4 

+ c|| v ^||^||v«||^||0||L»+q|v«|U4||v«|| ia ||0|U4 

<C\\V 2 9\\ L 2\\9\\ L 4 + C\\9\\ L 4 + C\\^p9\\ L 49\\ L 4 
+ C\\^\ L 4V^ L 49\\ L ~+C\\V^ L2 \\9\\ Li . 

It follows from (|3.42[) and Lemma 13.51 that 

||V 2 0|| L2 < C\\^9\\ L 2 + C\\V9\\ L 2 + C\\Vu\\ 2 L < + C 
<C\\^- P 9\\ L 2+C. 



9 pdx< 



dx 



+ 



p9dx 



<C\\V9\\ L 2+C\\^p9\\ L 2, 
which together with Poincare's inequality leads to 

\\9\\ L ,<C\\V9\\ L 2+C 

<C\\V9\\ L 2+C\\^p9\\ L2 . 

And can be estimated as follows, 

\\9\\ L ~<C\\V 2 9\\ L2 +C\\9\\ L 2 
<C\\^p9\\ L2 +C. 

Substituting (l335]) - (l3T58l ) to (13341) . we arrive at 

2jJ , „llv7flll2 



dx) + k||V0|| 2 2 



(3.54) 



(3.55) 



For the estimate for ||^||l4, we will follow the method used in Lemma 13.41 Let 9 
mfOdx, 



(3.56) 



(3.57) 



(3.58) 



Cv (J P\0f 

<C J p|^| 2 dx + C||^||L 2 ||V0|| L2 +C||V^|| i2 +C||^||L 2 (3-59) 
+ C||Vn|| L2 J p\9) 2 dx + C\\Vu\\ L 4^p9\\ L 2+C\\Vu\\ 2 L2 + j\\V9\\ 2 L2 , 

which together with the Gronwall's inequality completes the proof for Lemma 13.71 
As a corollary, the bounds for ||#||jy 2 and can be derived. 
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Corollary 3.8 Under the condition (13. ljl . it holds that forO<T< T* , 

sup {\\e\\ H 2 + ||0||£oo) < C. (3.60) 

0<t<T 

Proof. First, it follows from (13. 42ft . Corollary 13.61 and Lemma 13.71 that 

||V 2 #|| L 2 < C, (3.61) 

which implies that 

< C\\9\\ H 2 < C, (3.62) 
and hence completes the proof for Corollary! 



Up to now, we have get the bounds for ||p||l°° and ||0||^oo, which imply other nec- 
essary high-order estimates for the extension of the strong solution, according to the 
theorem proved in [6|. We sketch the proof for completeness. It is the same as that 
in [6] and the main idea follows that in |30j. 

Corollary 3.9 Under the condition (|3.1|) . it holds that for < T < T* 

fT 

sup \\w\\ H 2 + / \\(\\X7 2 w\\ 2 LP + \\Vw\\ 2 Lac )dt<C, pe(2,oo). (3.63) 

0<t<T JO 

Proof. By virtue of Proposition 12.11 and Lemma 13.51 

\\w\\ H 2 < C\\pu\\ L 2 < C, 
and by Sobolev embedding inequality, 

||Vio||l°o < C||Vio|| w i,p < C||pit||LP < C||?i||#i„ 
which implies (|3.63p . 

The next lemma is used to bound the density gradient and ||u||^2. 
Lemma 3.10 Under the condition (|3.1j) . it holds that 



sup (||p||wi,« + INI^a) < C. (3.64) 

0<t<T 



Proof. For 2 < p < q, |Vp| p satisfies the following equation 

(\Vp\P) t + div(\Vp\ p u) + (p - l)\Vp\ p divu 

+p\Vp\ p - 2 {Vp) tr Vu(Vp) + pp\Vp\ p - 2 Vp • Vdivu = 0. 



Hence, 
d 



^||Vp||iP < C(l + ||Vu|U«.)||Vp||iP + C||V 2 n|| LP 
dt 

< C (1 + ||V«||l- + ||Vw||oo) ||Vp||iP + C (||V 2 i.|| L p + ||V 2 u;|| LP ) ( 3 - 65 ) 

< C (1 + \\Vv\\ L oo + || VHloo) II Vp||iP + C||V 2 ^|| LP + c, 
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where we have used the fact 



\V 2 v\\ L v < C\\V(p8)\\ LP 



< C\\Vp\\lp\\0\\l°° +C\\V6\\ L p\\p\\l°° (3.66) 

< c\\v P \\ LP + c. 



To bound ||Vv||z,°o, we make use of the endpoint case of Proposition [2J] Lemma 
and (pT66]h 

HV^Hioo < C (l + ||Vw|| BMO (n) hi(e + UVulltyi.p)) 

< C (1 + + H^IM ln(e + \\Vv\\ w i, P )) (3.67) 

<C7(l + ln(e+||VH|Lp)). 
Substituting f|3.67j) into (|3.65p . we get that 

| + (3.68) 
<C(1 + ||Vw|U») ||Vp|| t » + Chi(e + ||Vp|| t »)||Vp||i» + C||V 2 u>||„ + C, 

which together with Gronwall's inequality gives that 

sup ||Vp|| L p < C. (3.69) 

0<T<T* 

Let p = q, then we get the bound of ||p||^i, g . 

Moreover, Let p = 2 in (|3,69p . then by Corollary 13.91 and (|3.66j) . 

||V 2 n|| L 2 < ||V 2 w|| L 2 + ||V 2 u|| L 2 < C + C||Vp|| z2 < C. (3.70) 

It completes the proof of Lemma 13.101 

In view of Lemma l3.U Lemma 13.101 it is enough to extend the strong solutions of 
(p,u,9) beyond t > T* . In fact, note that the generic constants C in Lemmas 13.11 - 
13.101 remains uniformly bounded for all T < T*, so the functions (p,u,6)(x,T*) = 
limt-^T* (p, u, 0)(x, t) satisfy the conditions imposed on the initial data at the time 
t = T*. Furthermore, standard arguments yield that pu, p6 G C([0, T];L 2 ), which 
implies 

(pu,p6)(x,T*) = lim (pii,p6) € L 2 . 

t— >T* 



Hence, 



-pAu- (p + X)Vdwu + RV( P e)\t=T* = y/p(x,T*) gi (x), 
kA8 + ||V« + (Vu) ir | 2 + A(divu) 2 | 4=T * = y/p(x, T*)g 2 (x), 



with 



p-y 2 (x,T*)(pu)(x,T*), for x E {x\p(x,T*) > 0}, 
0, for x £ {x\p(x,T*) = 0}, 



and 



92 (x) 



' p" 1 / 2 (x,T*)[p0 + Rp6dwu}(x,T*), for x G {x\p(x,T*) > 0}, 
0, for x £ {x\p(x,T*) = 0}, 
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satisfying gi,g2 G L 2 due to Lemma 13,101 and Corollary 13.81 Thus (p,u,9)(x,T*) 
satisfies (jl.lip and (| 1 . 12[) also. Therefore, one can take (p,u,9)(x,T*) as the initial 
data and apply Proposition 11.11 to extend the local strong solution beyond T* . This 
contradicts the assumption on T* . We thus finish the proof of Theorem 11.21 
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